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OPTIMIZED DEPLOYMENT OF PARTS IN A DISTRIBUTION NETWORK 

RELATED APPLICATIONS 

This application claims benefit under 35 U.S.C. § 119(e) of U.S. Provisional 
Application Serial No. 60/243,659, filed October 26, 2000, entitled "SYSTEM AND 
METHOD FOR OPTIMIZED DEPLOYMENT OF INVENTORY, OR RE- 
DISTRIBUTION OF EXISTING INVENTORY, ACROSS A MULTI-ECHELON 
DISTRIBUTION NETWORK." 

This application is related to U.S. Patent Application Serial No. 09/ , , 

entitled "REDISTRIBUTION OF PARTS IN A DISTRIBUTION NETWORK," by 
Rosa H. Birjandi, et al., Attorney's Docket No. 020431.0780, filed concurrently with 
the present application. 



TECHNICAL FIELD OF THE INVENTION 

This invention relates generally to the field of inventory distribution networks 
and more specifically to optimized deployment of parts in a distribution network. 
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BACKGROUND OF THE INVENTION 

Distribution networks may include one or more locations that receive parts 
from a vendor and distribute the parts within the distribution network in order to 
provide a customer with a product. The parts may be, for example, manufactured into 
a product within the distribution network. Distribution networks may include 
locations that both supply parts to and receive parts from other locations. 
Performance at each location is thus affected by the performance at its suppliers. As a 
result, maintaining an optimal inventory of parts at each location that best serves the 
customer while minimizing inventory costs poses a challenge for inventory managers. 
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SUMMARY OF THE INVENTION 

In accordance with the present invention, disadvantages and problems 
associated with inventory deployment and redistribution techniques are reduced or 
eliminated. 

According to one example of the present invention, a method for deploying 
parts is disclosed. Locations that include supply locations and demand locations are 
defined. A supply location supplies parts to a demand location. A demand is 
computed for each part at each location. An availability lead-time is estimated for 
each part at each location. A lead-time demand is computed for each part at each 
location using the availability lead-times for the part. A stock level is computed for 
each part at each location. A completely filled demand is determined from the lead- 
time demands and the stock levels, and a partially filled demand is determined from 
the lead-time demands and the stock levels. A coverage function for the parts at the 
locations is generated from the completely filled demand and the partially filled 
demand. 

Certain examples of the invention may provide one or more technical 
advantages. The present invention may be used to determine an optimized inventory 
deployment plan that describes the inventory at each location of a distribution 
network. The inventory deployment plan may optimize the ability of the distribution 
network to satisfy customer demand while conforming to business constraints. The 
inventory deployment plan may maximize the ability of the distribution network to fill 
orders, which may be calculated by minimizing the expected backorder of the 
distribution network. The present invention may be used to formulate a coverage 
function that is optimized to determine an optimized inventory deployment plan. 
Coverage may be used as a measure of customer service that describes the expected 
ability of each location to completely or partially fill a demand for a part. The present 
invention may be used to compute the expected number of partially and completely 
backordered demand for a part. 

The present invention may be used to calculate a net demand for a part at a 
location that accounts for dependent demands and independent demands. A 
dependent demand at a location describes the parts that the location supplies to other 
locations in the distribution network, and an independent demand at a location 
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describes the parts used at the location. Incorporating the independent and dependent 
demand into the demand may provide for a more accurate calculation of the demand. 
The present invention may be used to calculate a demand for a part at a location that 
takes into account the probability that the part is repaired and placed back into the 
inventory at the location. By taking into account the repaired parts, the calculation of 
the demand may be more accurate. 

The present invention may be used to calculate the availability lead-time for a 
part at any number of supply locations. The demand location may order a certain 
proportion of parts from the supply locations in a particular order. The computation 
of the availability lead-time takes into account the probability that a supply location 
supplies the part, given that no other supply location has supplied the part, which may 
provide a more realistic calculation of availability lead-time. The replenishment lead- 
time at a demand end point may be computed as the availability lead-time at its 
supplier plus the transfer lead-time from the supplier to the demand end point. 

The present invention may be used to calculate the expected number of 
backordered demand for a part at a location from the partially backordered and 
completely backordered demand. An equivalence relation between maximizing the 
coverage function and minimizing the sum of backorders may be determined. 

Other technical advantages may be readily apparent to one skilled in the art 
from the figures, descriptions and claims included herein. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of the present invention and its features 
and advantages, reference is now made to the following description, taken in 
conjunction with the accompanying drawings, in which: 

FIGURE 1 illustrates an example distribution network for deploying and 
redistributing inventory of one or more parts among one or more locations; 

FIGURE 2 illustrates an example system that generates optimized inventory 
deployment and redistribution plans; 

FIGURE 3 illustrates an example method for deploying and redistributing 
inventory of one or more parts among one or more locations; 

FIGURE 4 illustrates an example method for calculating a demand for one or 
more parts at one or more locations; 

FIGURE 5 illustrates an example method for estimating the availability of one 
or more parts at one or more locations; and 

FIGURE 6 illustrates an example method for generating a coverage function 
for one or more parts at one or more locations. 
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DETAILED DESCRIPTION OF THE DRAWINGS 

FIGURE 1 illustrates an example distribution network 20 for deploying and 
redistributing inventory of one or more parts among one or more locations 22. 
Distribution network 20 includes locations 22 that distribute parts throughout 
distribution network 20. A part may comprise, for example, a product, a portion of a 
product, a device used to manufacture a product, or any other suitable item that may 
be distributed from one location 22 to another location 22 in distribution network 20. 

In one embodiment, locations 22 include a central location 22a and one or 
more warehouse locations 22b-d. Although central location 22a and warehouse 
locations 22b-d are illustrated, distribution network 20 may include any suitable 
number of central locations 22 and warehouse locations 22. Each location 22 may 
comprise a supply location and/or a demand location. A supply location supplies a 
part to a demand location, and may supply the part in response to an order for the part 
sent from the demand location. For example, warehouse location 22b supplies parts 
to warehouse location 22d, and warehouse locations 22b-c supply parts to location 
22d. A location 22 may comprise both a demand location and a supply location. For 
example, warehouse location 22b receives parts from central location 22a and 
supplies parts to warehouse location 22d. A supply endpoint such as central location 
22a receives parts from one or more external supplies 24, for example, a vendor, and 
distributes the parts to warehouse locations 22b-d. A demand endpoint such as 
warehouse location 22d provides parts to one or more external demands 32, for 
example, a customer. 

Warehouse locations 22b-d may include supply operations 26b-d and/or repair 
operations 28b-d. A supply operation 26 sends an order for a part to a supply 
location, which in response sends the part to supply operation 26. A repair operation 
28 may receive a broken part from supply operation 26 and send the broken part to a 
repair center 30. Repair center 30 repairs the part and sends the repaired part to, for 
example, central location 22a or back to supply operation 26b. Alternatively, repair 
operation 28d may receive a broken part from supply operation 26d, repair the part, 
and send the repaired part back to supply operation 26d. 

The inventory for each part at each location 22 is monitored, continuously or 
periodically. In response to the inventory falling below a predetermined level, an 
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order is placed to bring the inventory position back up to a target level such as an 
optimized inventory level. A method for deploying and redistributing inventory of 
one or more parts among one or more locations to achieve optimized inventory levels 
is described in more detail with reference to FIGURE 2. 

FIGURE 2 illustrates an example system 34 that generates optimized 
inventory deployment and redistribution plans. An inventory deployment plan 
describes a distribution of parts among locations 22 of distribution network 20, and an 
inventory redistribution plan describes a manner of redistributing parts to satisfy an 
inventory deployment plan. Deployment may occur independently of redistribution. 
That is, inventory may be deployed, without ever being redistributed. Additionally, 
redistribution may redistribute parts according to any suitable inventory plan. 

System 34 may include a computer system 35, a server 36, and a database 37, 
which may share data storage, communications, or other resources according to 
particular needs. Computer system 35 may include appropriate input devices, output 
devices, mass storage media, processors, memory, or other components for receiving, 
processing, storing, and communicating information according to the operation of 
system 34. As used in this document, the term "computer" is intended to encompass a 
personal computer, workstation, network computer, wireless data port, wireless 
telephone, personal digital assistant, one or more microprocessors within these or 
other devices, or any other suitable processing device. 

Server 36 manages applications that generate optimized inventory deployment 
and redistribution plans. Server 36 includes one or more software components such 
as a pre-processing module 38 and a solver 39. Pre-processing module 38 manages 
input and output operations, and computes a net demand and a replenishment lead- 
time for each part at each location 22. Pre-processing module 38 computes a demand 
over lead-time, or lead-time demand, related to a number of parts in a pipeline. Pre- 
processing module 38 also generates mathematical formulations, which are 
transmitted to solver 39 for solving. 

Pre-processing module 38 may include a deployment module 40 and a 
redistribution module. Deployment module 40 may be used to generate a coverage 
function that describes the distribution of parts among locations 22. Solver 39 
optimizes the coverage function to determine an optimized distribution of parts. 
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Solver 39 may comprise any suitable mathematical programming solver such as 
CPLEX by ILOG, INC. Redistribution module may be used to generate a transfer 
function that describes the cost of transferring parts among locations 22. Solver 39 
optimizes the transfer function to determine an optimized manner of redistributing 
parts. As noted above, deployment may occur independently of redistribution. That 
is, inventory may be deployed, without ever being redistributed. Additionally, 
redistribution may redistribute parts according to an inventory plan generated by 
deployment module 40 or according to any suitable inventory plan. 

Database 40 stores data that may be used by server 36. Data may include, for 
example, the history of the demand for each part at each location 22, the lead-time 
required to transport a part from one location 22 to another location 22, and the 
maximum space capacity at location 22. Computing system 35 and database 40 may 
be coupled to server 36 using one or more local area networks (LANs), metropolitan 
area networks (MANs), wide area networks (WANs), a global computer network such 
as the Internet, or any other appropriate wired, optical, wireless, or other links. 

FIGURE 3 illustrates an example method for deploying and redistributing 
inventory of one or more parts among one or more locations 22. Deployment may 
occur independently of redistribution. That is, inventory may be deployed, without 
ever being redistributed. Additionally, redistribution may redistribute parts according 
to the inventory deployment plan described in connection with FIGURE 3, or 
according to any suitable inventory plan. 

Processing module 38 initiates the method at step 46 by defining a number 1, 
2, . . ., i 9 . . I of parts and a number 1,2,.. . . J of locations 22. For example, 
j=l, 2, 3 and 4 and refer to warehouse locations 22a-d, respectively. At step 48, data 
is accessed from database 37. Data may include, for example, a demand history of 
each part at each location 22. The demand history may describe the parts that each 
location 22 requires. Data may include repair history that may describe the capability 
of each location 22 to repair a part. Data may include the paths that may be used to 
transfer parts between locations 22, along with the costs associated with transporting 
parts along the paths. Data may include the cost of purchasing a part, the cost of 
storing a part in the location as a percentage of the purchase cost for the part, and a 
cost associated with ordering a part. 



ATTORNEY' 
020431.0787 




:ket 



# 



'ENT application 



At step 50, a demand for each part at each location 22 is calculated. The 
demand may include a dependent demand and an independent demand. A dependent 
demand at location 22 describes the parts that location 22 supplies to other locations 
22. An independent demand at location 22 describes parts used at location 22. The 
demand at location 22 may account for the probability that a part is repaired and 
placed back into the inventory at location 22. Demand may be calculated by starting 
at a demand endpoint and ending at a supply endpoint of distribution network 20. A 
method for calculating a demand for a part at each location 22 is described in more 
detail with reference to FIGURE 4. 

A replenishment lead-time for each part at each location 22 is calculated at 
step 52. The replenishment lead-time for a part at location 22 describes the time 
required for location 22 to receive the part from another location 22. The 
replenishment lead-time may be computed by starting at a supply endpoint and ending 
at a demand endpoint. An availability lead-time for each part at each location 22 is 
estimated at step 54. The availability lead-time at a location 22 describes a waiting 
time due to back order at location 22 plus the transfer lead-time from a supplier to 
location 22 and the replenishment lead-time for the supplier of location 22. A method 
for estimating the availability lead-time of a part at location 22 is described in more 
detail with reference to FIGURE 5. 

A coverage function is determined at step 56. The coverage function 
describes the expected ability of a location 22 to completely or partially fill an order 
for a part, and may be determined from the demand, availability lead-time of the part, 
and inventory level for the part at location 22. The coverage function may be 
described using the expected backorder of the part at location 22. A method for 
determining a coverage function is described in more detail with reference to 
FIGURE 6. Solver 39 optimizes the coverage function at step 58. Optimizing the 
coverage function may be accomplished by minimizing the expected backorder. At 
step 60, an optimized inventory for each part at each location 22 is determined from 
the optimized coverage function. At step 61, solver 39 reports the optimized 
inventory for each part at each location 22. As noted above, deployment may occur 
independently of redistribution. That is, inventory may be deployed, without ever 
being redistributed. 
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At step 62, redistribution module determines whether redistribution is required 
by calculating an excess and deficit for each part at each location 22 from the actual 
inventory and the optimal deployment. Redistribution of the inventory may be 
required if, for example, the actual inventory at each location 22 does not match the 
optimized inventory calculated for each location 22. As noted above, redistribution 
may redistribute parts according to the optimized inventory reported at step 61, or 
according to any suitable inventory plan. If redistribution is not required, deployment 
module 40 proceeds to step 63 to report any excess inventory and a recommendation 
to not perform a redistribution of parts. After reporting the result, the method is 
terminated. If redistribution is required, redistribution module proceeds to step 64 to 
check the transitions between locations 22. The transitions describe paths that may be 
used to transfer parts from one location 22 to another location 22. 

At step 66, a transfer function describing the transfer of parts between 
locations 22 is optimized. Minimizing the total costs associated with transporting the 
parts may optimize the transfer function. A method for determining optimized 
transfer plans for transferring parts between locations 22 is described in more detail 
with reference to FIGURE 7. At step 70, the optimized transfer plans, the resulting 
inventory levels, and possible excess in inventory of parts in the network are reported. 
After reporting the result, the method is terminated. 

FIGURE 4 illustrates an example method for calculating a demand for one or 
more parts at one or more locations 22. Deployment module 40 initiates the method 
at step 80 by selecting a part i. A location j is selected at step 82. Location j may be 
selected such that the demand at a demand endpoint is calculated first, and the 
demand at a supply endpoint is calculated last. 

At step 84, an independent and a dependent demand for part i at location j is 
determined. The independent demand for part i at location j may be represented by 
Vij. The dependent demand for part i at location j may be represented A,,*, where A: is a 
demand end point for location j. At step 86, the repair capability ry for part i at 
location j is determined. The repair capability ry may be determined from the 
proportion of demand for part i at location j that is repairable at location j. Demand is 
calculated at step 88. Starting with demand end points j 9 demand Xy for part / is equal 
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to its independent demand. For any location j that is not a demand end point, the 
demand for part i may be calculated using Equation (1): 



At step 92, deployment module 40 determines whether there is a next part for 
which a demand is to be determined. If there is a next part, deployment module 40 
returns to step 80 to select the next part. If there is no next part, deployment module 
40 proceeds to step 94 to output the calculated demand for each part at each location. 
After outputting the demand, the method is terminated. 

FIGURE 5 illustrates an example method for estimating the availability lead- 
time of one or more parts at one or more locations 22. Deployment module 40 
initiates the method at step 102 by selecting a part i. A demand location j is selected 
at step 104, and a supply location 4 is selected at step 106. The supply location 4 
may be selected from a prioritized list of n of supply locations //, /„. For each 
supply location 4, the list may describe a proportion Cukj of a demand for part i at 
demand location j that is scheduled to be satisfied by supply location 4, a probability 
amy that part i is filled at supply location 4 for demand location j 9 and a lead-time T ak j 
for a part i to flow from supply location 4 to demand location j. Demand location j 
and supply location 4 may be selected such that the availability lead-time at a supply 
endpoint is calculated first, and the availability lead-time at a demand endpoint is 
calculated last. 

At step 108, a probability Pu k j of a supply location 4 filling an order for part i 
placed by a demand location j, given that the order is not filled by another supply 
location, is calculated. The probability Paij for supply location // may be computed 
using Equation (2): 

At step 110, deployment module 40 determines whether there is a next supply 
location l k . If there is a next supply location, deployment module 40 returns to step 
106 to select the next supply location. The probability Pay of the next supply location 
4 filling an order for part i placed by demand location j 9 given that the order is not 
filled by another supply location, may be computed at step 108 using the process 
described by recursive Equations (3): 




'ik 



(i) 



k is a demand point for j 
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p uj = a u k j C i k j 
where 

C =C (3) 

Qy=C, 7ty+ (l-a, w )c, w fork>l 

If there is no next supply location at step 110, deployment module 40 proceeds to step 
112 to output the probabilities of the supply locations 4 fulfilling an order placed by 
demand location j. 

5 At step 113, an availability lead-time 7^ for each location j is calculated. 

H* Availability lead-time Ty may be calculated according to recursive Equation (4): 



ru m 



k i'k 



P n k J (4) 



The expected number of completely backordered demand B c may be described by 
Equation (5): 

5 10 £»O^^JxP(jf|jiJ (5) 

m 

US / I \ 

p where P[X\/i 9 )=e ^ii^* / x\ is the Poisson probability mass function for the 

distribution of demand with mean ^ ja,* represents the mean number of parts i in the 
pipeline at supply location k. The expected number of partially backordered demand 
B p for part i at location j may be described by Equation (6): 

1 5 EBO p =(l^ z )f j ( x - SiJ )p(x\ ix , ) (6) 

x=S u 

The expected number of backorders EBO (Sij) having the stock level Sy of part / at 
location j may be defined using Equation (7): 

EBO(Sij) = Z ± xP(x\ iiiJ ) H l^ z )f d ( x -S ij )p(x\ n , ) 

X=Sy X= Sjj 

= Z(x-*S ij )p(x\ii u ). (7) 



x=S n 



20 At step 114, the replenishment lead-time for part i at demand location j is 

calculated. The replenishment lead-time S,y for part i at location y may be calculated 
using Equation (8): 
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(8) 



where xg represents the repair lead-time for part i at demand location j. The lead-time 
demand \ig of a part i at demand location j is estimated at step 116. The lead-time 
demand is related to a number of parts in a pipeline. The lead-time demand may be 
estimated using Equation (9): 



At step 118, deployment module 40 determines whether there is a next 
demand location. If there is a next demand location, deployment module 40 returns to 
step 104 to select the next demand location. If there is no next demand location, 
deployment module 40 proceeds to step 120 to determine whether there is a next part. 
If there is a next part, deployment module 40 returns to step 102 to select the next 
part. If there is no next part, deployment module 40 proceeds to step 122 to report the 
lead-time demand of each part at each location 22. After reporting the lead-time 
demand, the method is terminated. 

FIGURE 6 illustrates an example method for generating a coverage function 
for one or more parts at one or more locations 22. Deployment module 40 initiates 
the method at step 132 by selecting a location j. A part i is selected at step 134. 

At step 136, a completely filled demand D c for part i at location j is calculated 
at step 136. A completely filled demand D c may be described by Equation (10): 



where P[X\/J 0 .)=e ^ J \x 0 x / x\ is the Poisson probability mass function for the 

distribution of demand with mean yaj. A partially filled demand D p for part i at 
location j is calculated at step 138. The partially filled D p demand may be described 
by Equation (11): 



where % is the percentage of partial fill allowed for a part. At step 140, a coverage 
function for part i at location j is determined. The coverage function for part / at 
location j describes the expected proportion filled demand for part i of location j\ and 
may be expressed using Equation (12): 



(9) 




(10) 



jc=0 




(11) 
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\txP(x\ Mu )+ Z ±(s y -iy{x\A/^ = (12) 
&(x- Z(S y - l))P(X\ju ij ) + % (S U -1)}/M, 

x=0 

At step 142, deployment module 40 determines whether there is a next part. If 
there is a next part, deployment module 40 returns to step 134 to select the next part. 
If there is no next part, deployment module 40 proceeds to step 144 to determine 
whether there is a next location. If there is a next location, deployment module 40 
returns to step 132 to select the next location. If there is no next location, deployment 
module 40 proceeds to step 146 to determine the coverage function for the number of 
parts at the number of locations. The coverage function may be expressed as the 
weighted average of coverage for the parts at locations. The coverage function for the 
parts at the locations may be expressed by Expression (13): 

t It A £ l* - xs v M4*, )+ X s v ) / ± n, (13) 

y = l [ , = 1 X =0 J 1=1 

where (3, represents a weight of part i, which may be based on an importance measure 
of part i. At step 148, constraints for the coverage function may be defined. 
Constraints may include, for example, the following: 

a. The weighted average of coverage for the parts at each location j is 
greater than or equal to the coverage target co, at location j 9 which may be expressed 
by Expression (13a): 

>co,, V/ (13a) 

1=1 

b. The coverage for each part i at each location j is greater than or equal 
to the coverage target for part i at location j 9 which may be expressed by Expression 
(13b): 

s 

Y,[x-zS 8 ]p(x\ Mu )+zSy >a y H,,V/,V/ (13b) 

JC=0 

c. The number of new purchases for a part i, Xy, at location j, which may 
be expressed by Expression (13c): 



1=1 \x=0 
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^-fo+r^-itf, vi.vy (13c) 

where y y represents a proportion of a demand for a failed part i , and (x) + = max (0^). 

d. The inventory investment at each location j is less than or equal to the 
inventory investment limit Invj for location j, which may be expressed by Expression 
13d): 

j^cfaj+Y^j-Y^InVj, V/ (13d) 
1=1 

where Q represents a purchase price for part i, and Yy represents an on-hand inventory 
for part i at location j. 

e. The overall inventory investment is less than or equal to an overall 
inventory investment limit Inv, which may be expressed by Expression (13e): 

'£c g X 9 *Inv (13e) 
/=/ 

f. The overall inventory cost is less than or equal to overall budget B, 
which may be expressed by Expression (13f): 

iltskMj ~M, + J3BO, + £ ±C iXy <B (13f) 

where hy represents a holding cost per unit of part i at location j 9 and ki 
represents an order cost for part i. 

g. The cost at each location j is less than or equal to a budget Bj at 
location j 9 which may be expressed by Expression (13g): 

2^,(5, -Mi+EBOJ + kfa+CiSi+rfa -Y 9 Yl*Bj, Vj (13g) 

1=1 

h. The total volume occupied by parts at each location j is less than or 
equal to the volume capacity limit Vj at location j 9 which may be expressed by 
Expression (13h): 

i>,.s,<r,., vy (i3h) 

i. The stock levels are integers, which may be expressed by 
Expression (13i): 

Sy are integers ( 1 3 i) 
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At step 150, the coverage function is converted to a backorder function that 
corresponds to expected backorders, and the constraints are expressed in terms of 
backorders. Using the backorder function may provide for a simpler optimization 
process. The backorder function may be Expression (14): 

±±EBO y (5,) (14) 
j=\ i=i 

Minimizing the backorder function is equivalent to maximizing the coverage function. 
The constraints may be expressed by Expressions (14a): 

S V >S%, Vi,V/ (14a) 

Y J EBO ij {S ij )<(X-m j )^y, vy 

it&yiSy - MiJ +EBO y ) + k i X iJ ]} + f i f j C i x iJ <B 

j=\ 1=1 7=1 1=1 

£ [hy (Sy - My + EBOy ) + kfa + C, (s, + - 7 # )\ < Bj , Vj 
1=1 

•tv,.5,<^, Vy 

1=1 

5^- are integers 

where represents the minimum stock level that guarantees the minimum target 

coverage. The high degree of non-linearity of the constraints may be reduced by 
replacing the minimum target coverage constraint for each part at each location 
expressed by Expression (14a) with an equivalent constraint expressed as 
S y > Sy , Vz, V/. The definition of the number of new purchases expressed by 
Expression (14b) may be replaced with a relaxed constraint expressed as 

Xy>S u+ \ ri Xy-YylVi. 

Maximization of the coverage function may be shown to be equivalent to 
minimizing the backorder function. The weighted average of coverage for 
distribution network 20 may be expressed by Expression (14b): 

Hfifi 9 (S 9 ) (14b) 

y=l i=l 

where QijiSij) represents the coverage for part i at location j defined by Equation (14c): 
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0y(S ij ) = (£ j xP(X\ Mij ) + X ^S 0 P(X\m„)}/h„ (14c) 

The backorder function describes the total expected number of backorders for 
distribution network 20 and may be expressed by Expression (14d): 

HeB0 9 (S 9 ) (14d) 
y=i i=i 

The following proposition may be established: 

Proposition 1: Vector Sj=(Sjj, ...,Sjj ) satisfies the performance constraint expressed 
by Expression (14e): 

X//^(5,)>^in, (14e) 

/=1 i=l 

if and only if vector Sj satisfies an expected backorders constraint expressed by 
Expression (14f): 

£eB0 9 {S v ) < (1 - a>j) £ V„ (14f) 

j'=l i=l 



A relationship between expected backorder EBO,y (Sy) and performance QyiSy) may be 
established according to Equations (14g): 

e v {S u ) = { MiJ - f J xP<<X\M i j) + X t d S iJ P{X\ f x u )}l V i iJ (14g) 
0 y (S y ) = {juy - X (x - zS 9 )P(X I juy )} / » v 

x=S u 

By (Sy ) = { M y - EBO, (S y )}///,= 1 - {£SO, (5, )} / ^ 

^o,(5,) = //,[i-e,(s,)] 

For any set of S/j that satisfies the performance constraint expressed by 
Expression (14f), the following Expressions (14h) may be shown: 

Y j EBOy(Sy)<(\-0> j )f jV iy 

i=\ i=l 

Z ft, [i - o § (Sy , (1 4h) 



ATTORNE^^bOCKET 

020431.0787 



ATENT APPLICATION 



18 



ii 11 

Z P(> <~ Z (Sy ) <Z Mij ~ Oj Z ^ 

i=\ 1=1 r=l /=1 

/ / 

-Z^W^-^Z^ 

1=1 »=1 
Z>"A(^)^;Z^ 

r=l 1=1 



The stock level Sy in question satisfies the location performance constraint expressed 
by Equations (14g). The steps may be reversed to prove the converse. 

The following proposition describing the relationship between performance 
and coverage may be established: 

Proposition: Maximizing system-wide coverage is equivalent to minimizing the total 
system-wide backorders. 

The proposition may be established according to Equations (14i): 
Maximize 

Z2>#°# (s 9 ) = ZZ 1 -^- (s„) (140 

j=\ i=i y=i i=i 

= 1 - Z Z EBO iJ ) = -tYZ EBOy (Sy )]-l 

y=l i=\ y=l i=l 

Maximize - £ £2?<9, (S, J>/ - / 
== Minimize ]T £ £50-. (S.. j 

7=1 i=l 

At step 154, the objective function that measures expected backorder as 
expressed by Expression (13) may be linearized. To linearize the objective function 
and constraints, the non-linear terms of the objective function and constraints may be 
approximated by linear terms. The non-linear terms are discrete and convex, so a 
first-order linear approximation using the finite difference for two neighboring 
discontinuous points may be used to approximate each non-linear term. Each non- 
linear term in the objective function and the constraints is replaced with a continuous 
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variable ty, and a linearization constraint that describes the under estimation at points 
of discontinuity is added to the constraints. 

The linearized objective function may be expressed by Expression (15): 

it*. < i5 > 

The linearization constraint may be expressed by Expression (16): 

't, > m, (x - X 9 ) + 6, , V5, < X, < S upper , ViJ (16) 

where m, = p(x > X^), b i} =p(x>X ( .\ i x 0 \x - X„) + EBO^X,. + l), and S upper is the 

upper bound on the inventory for part i at location j. Other constraints may be 
expressed by Expressions (16a): 

s y >s^, v/ = i,.../ 1 ,y 7 - 



i=\ x=0 



/ 2>, r >co,. , V/ 

X 9 +r i h-Y 0 , Vi = 7,..,7 (16a) 

1=1 

j=\ 1=1 

£c,.(s,+^-^)^"v,, V/- 

1=1 
1=1 

5 /y are integers, V/ = l,..,/ 2 , V/ 

After linearizing, deployment module 40 solves the resulting mixed integer- 
programming problem and determines the optimal stock levels for each part at each 
location. 
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An objective function measuring the system-wide total cost may also be 
defined at step 154. The total cost function may be expressed by Expression (16b): 

ZBW -M,+EBO,) + k,X, +C,x # ] (16b) 

7=" '=1 

The constraints for the total cost function may be expressed by Expressions (16c): 
S u >S%, Vi=l,A 

[£ {Z A - V>(X I /<„, ) + xS,j }]/][>,.> co , , Vy 

^ >5, + ~Y 0 } Vi = l,..,I (16c) 

Zv f 5, ZVj, Vy 

5" are integers, Vz = 1, . /, V/ 

The total cost function and the constraints may also be linearized at step 154 in 
order to allow the objective function to be optimized by solver 39. The total cost 
objective function, as expressed by Expression (16b), may be linearized according to 
Expression (17d): 



j i 



ZZtW -Mi + (16d) 

y=l ,=l 

The constraints may be linearized according to Expressions (16e): 

t t > m t) (S -X„) + b y , VS? < Xj < S upper , V/j 
Sy >S$, Vi = l,A/„Vy 

t^>(l-« y ;tu,. Vy >S u+ \ ri Xy-Yy} Vi = l,..,I 

1=1 1=1 

(16e) 

£v,.S, <F,, Vy 

1=1 

5// are integers, V/ = 1,..,/ 2 ,V/ 

After linearizing, deployment module 40 solves the resulting mixed integer- 
programming problem and determines the optimal stock levels for each part at each 
location. 

Although an example of the invention and its advantages are described in 
detail, a person skilled in the art could make various alterations, additions, and 
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omissions without departing from the spirit and scope of the present invention as 
defined by the appended claims. 



